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An analytic solution for Helfrich spontaneous curvature membrane model (H. Naito, M.Okuda 
and Ou-Yang Zhong-Can, Phys. Rev. E 48, 2304 (1993); 54, 2816 (1996)), which has a conspicuous 
feature of representing the circular biconcave shape, is studied. Results show that the solution 
in fact describes a family of shapes, which can be classified as: i) the flat plane (trivial case), 
ii) the sphere, iii) the prolate ellipsoid, iv) the capped cylinder, v) the oblate ellipsoid, vi) the 
circular biconcave shape, vii) the self-intersecting inverted circular biconcave shape, and viii) the 
self-intersecting nodoidlike cylinder. Among the closed shapes (ii)-(vii), a circular biconcave shape 
rf\ ' is the one with the minimum of local curvature energy. 
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I. INTRODUCTION 



Why red blood cells under normal physiological conditions take the circular biconcave shape (CBS) but not the 

""O ■ spherical shape aroused the long standing curiosity of human being since its first discovery in the seventieth century. 

Today, physicists ascribe it to the minimization of bending energy of flexible lipid bilayer membrane consisting of 

amphiphilic molecules. The first successful model revealing the morphology was due to Canham (fl in 1970. However, 

his theory suffered from the shortcoming that the membrane was assumed to consist of two identical labile surfaces, 

and both chemical and physical environments between two sides of the membrane were assumed to be identical too. 

To overcome the shortcoming, Helfrich in 1973 introduced a phenomenological parameter, namely the spontaneous 

qq . curvature, in describing the more realistic situations: the asymmetry of the two leaflets of the membrane and the 

£f) ' chemical or/and physical differences between the interior and exterior membrane B. As expected, this model gave 

more abundant shapes than that of Canham. Based on the numerical integration technique |3j, Helfrich spontaneous 

curvature model yielded a catalog of axisymmetric vesicle shapes: the CBS, the prolate and oblate ellipsoid, etc. More 

than ten years later, the general equilibrium equation was derived by performing the variation of the Helfrich energy 

qv ' functional Q . The first triumph of the equation was the prediction of the existence of Clifford torus shape membrane 

p] and its subsequent experimental verification J6|. In 1993, the general equation in axisymmetric case was obtained, 

and it is a complicated third order nonlinear differential equation |J7j . An analytic solution capable of representing 

the CBS H was immediately obtained. In the same year, the first integral of the third order nonlinear differential 

equation was found; then the equation reduced to be a second order one involving a constant of integration C |J . An 

interesting fact is that the analytical solution requires a nonvanishing value of C ||[| , whereas only the special case 

of the second order differential equation with C — has been well studied by numerical methods pQl. I n 1996, the 

solution was used to explain the experimentally observed polygonal shape transformation of the CBS |ll[| . Studies in 

this paper will show that this solution actually represents a family of shapes in this solution, but in the family, the 

• i— i . CBSs have the lower energies and a CBS has the lowest. For convenience, we will call this solution the CBS solution 

y\ ' hereafter. 

The existence of the analytical CBS solution in Helfrich model proves to be remarkable, because the analytical 
solutions to a nonlinear theory are rare and precious |13|. In very special cases of Co = dp = A = (cf. following 
Eq.(l)), the Helfrich energy functional reduces to that for constant curvature surfaces and Willmore surfaces p3[, 
and all known analytical solutions having physical applications in membrane shapes are Delaunay surfaces, sphere, 
torus J13J , and no more E3. Therefore a systematic study of the CBS solution is necessary. We will find that there 
are eight types of shapes involved in the solution and the enclosed shapes can be grouped into the prolate or oblate 
ellipsoid branches. These two branches are bifurcated from a sphere. 

The article is organized as follows. In section II, how to obtain the CBS solution from the Helfrich model is outlined. 
In section III, all typical shapes contained in the CBS solution are plotted and their parameterizations are presented. 
In section IV, a systematically analysis of the CBS solution is given. In the shape family, there is an shape with 
minimum energy and in section V this shape is found from the the scale invariance of the local curvature energy. In 



section VI, a comparison of our results with the previous experimental and theoretical results is given. In the final 
section VII, a brief conclusion is given. 

II. HELFRICH SPONTANEOUS CURVATURE MODEL AND ITS CBS SOLUTION 

Equilibrium shapes of phospholipid vesicles are assumed to correspond to the minimum of the elastic energy of the 
closed bilayer membrane. The energy functional of Helfrich spontaneous curvature model reads as B 

F = ~/c / (ci + c 2 - c ) 2 dA + Sp J dV + A I dA, (1) 

where dA and dV are the surface area and the volume element for the vesicle, respectively, k is an elastic modulus, 
c\ and C2 are the two principal curvatures of the surface, and cq is the spontaneous curvature to describe the possible 
asymmetry of the bilayer membrane. When cq is zero Helfrich model reduces to Canham model H. The Lagrange 
multipliers Sp and A take account of the constraints of constant volume and area, which can be physically understood 
as the osmotic pressure between the ambient and the internal environments, and the surface tension, respectively. 
The general equilibrium shape equation is B 

Sp - 2XH + k{2H + c a ){2H 2 - 2K - c H) + 2k\j 2 H = 0, (2) 

where SJ 2 = -j=di{g l: ' y/gdj) is the Laplace-Beltrami operator, g is the determinant of the metric gij and g^ = {gij)~ l , 
K — c\Ci is the Gaussian curvature and H = — (l/2)(ci + c^) is the mean curvature. Assuming that the shape has 
axisymmetry, the general shape equation Eq.(2) becomes a third order nonlinear differential equation [Q 
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where r is the distance from the symmetric z axis of rotation, V'( r ) is the angle made by the surface tangent and the 
r axis as shown in FIG.l. The positive direction of the angle is that the angle measured clockwise from r axis. It 
is contrary to the usual mathematical convention; therefore, the mean curvature H is —\/2{su\ip /r + dsimp/dr), in 
which c\ — smijj/r denotes the principal curvature along the parallels of latitude, and C2 = d simp /dr denotes that 
along those of meridian. It is worthy to mention that the spontaneous curvature cq carries a sign. When the normal 
of a surface change its direction, c\, ci and cq must change their signs simultaneously. Keeping the directions of r and 
z as the usual, we have consequently 

dz/dr — — t&iiijj(r) 

z(r) - z(Q) = - j* t&nip(r)dr (4) 

n = (sini/; cos</>, sin-0 sin0, cost/'), 

where n denotes the normal of the surface and 4> the azimuthal angle. 

For self-consistence, the positive direction of arclength s along the contour in r — z plane is necessary to start from 
the north pole of the shape. If so, we see that the parameterization for sphere is sin-0(r) = r/Ro where Rq > is the 
radius of the sphere, and our sign convention used in this paper is then compatible with that used in most previous 
works pMGll . The third order nonlinear differential equation (3) can be simplified to be a second order one M 
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where C is a constant of integration. It is still a rather complicated equation and does not belong to any known type 
of well-studied differential equation in mathematics. 



Under the conditions that both the surface tension A and the osmotic pressure difference Sp are zero, i.e., 

A = Sp = 0, (6) 

an analytic CBS solution for the Eq.(3) [§0, or Eq. (5) with C = 2c , is 

sin-0 = r/Ro + cor lnr, or equivalently, sin?/' = cor ln(r/r m ), (7) 

where i?o is an arbitrary constant and r m = exp(— l/(co-Ro))- When we first obtained this solution pjll[, we 
concentrated on the fact that it can be used to represent the CBS of usual red blood cells pj| . In fact, by adjusting 
parameter cq in the interval (— oo,oo), the formula (ufi can give a family of shapes, which is what we are going to 
analyze in detail. For sake of convenience, we mainly use the first form of the CBS solution (pi) with Rq = 1. 

III. ALL POSSIBLE SHAPES IN THE CBS SOLUTION 

Before presenting the shapes, we would like to make two comments on how to characterize the shapes. First, for 
each shape, we will give both the parameterization and the interval of r e [0, oo) in which the shape appears. This is 
because that there may be different shapes represented by the same parameterization in distinct intervals of r e [0, oo) 
as long as in the intervals | smip(r)\ < 1 are satisfied. Second, even the scale invariance will be fully studied in section 
V, we will follow the usual usage |10| to give for each shape the scale invariant cor s , in which r s = ^JA/An with 
A denoting the area of the surface of the shape. Unless specially mentioned and discussed, that the normals of the 
surface point outwards for closed shapes is always implied. 

All possible types of shapes represented by Eq. (Q) are listed in the following. 

i) The flat plane (trivial case) with cor s = 0: 

sinV> = 0, re [0,oo]. (8) 

ii) The sphere with cor s = 0: 

simP = r/Ro, re [0,i? ]. (9) 

hi) The prolate ellipsoid. A typical shape with c$r s = —0.72 is shown in FIG. 2(a) and one of its parameterization 

is 

sin V = r - 0.6r In r, r e [0, 1] . (10) 

Another parameterization representing the same shape is sin?/; = r — 1.85rlnr, r E [0,0.324]. The parameterization 
sin-0 = r + 4.063r lnr, r £ [0, 0.148] represents the same shape but with inwards pointing normal. 

iv) The capped cylinder. A typical shape with cor s — —2.06 is shown in FIG. 2(b) and one of its parameterization 

is 

sinV> = r-0.99rlnr-, re [0,1]. (11) 

Another parameterization representing the same shape is sin?/; = r — l.Olrlnr, r e [0,0.980]. The parameterization 
sin?/; = r + 3.5913r lnr, r e [0, 0.275] represents the same shape but with inwards pointing normal. 

v) The oblate ellipsoid. A typical shape with cor s = 0.46 is shown in FIG. 2(c) and its parameterization is 

sin'0 = r + 0.5rlnr, re [0,1]. (12) 

vi) The CBS. A typical shape with cor s = 1.51 is shown in FIG. 2(d) and its parameterization is 

s'mip = r + 1.8rlnr, re [0,1]. (13) 

vii) The self-intersecting inverted CBS. A typical shape with c$r s = 2.72 is shown in FIG. 2(e) and its parameteri- 
zation is 

sin-0 = r + 3.2rlnr, re [0,1]. (14) 



To note that in this case, the outside of this shape can be defined but we can not let the normal always point outwards 



using a single parameterization (14). The normal of the parameterization ( |14| ) points outwards only at the surface of 
the torus part. 

viii) The self-intersecting nodoidlike cylinder. This situation was discussed in previous paper p2| . We add a typical 
figure, FIG. 2(f), in this paper for completeness. The cylinder is infinitely long with periodic packing, along the 
rotational axis, of a basic unit in which self-intersecting occurs once. We plot the basic unit only. We use the surface 
of the basic unit to calculate r s and the shape has cqt s — 3.28. Its parameterization is 

sin7/> = r + 3.60rlnr, re [0.301,1]. (15) 

The normal of this parameterization points outwards only at the surface of the torus part also. To note that the same 
parameterization simp = r + 3.60rlnr in interval r G [0, 0.257] gives a capped cylinder but its normal points inwards. 
These eight shapes consist of all possible types of shapes containing in the solution. From the parameterizations 
of the these shapes, one may have noticed two facts that the spontaneous curvatures of all shapes are within a very 
narrow domain instead of infinite domain of cq S (— oo, oo), and a single shape may have different parameterizations. 
Explanations of these facts will be given in next section by both qualitative and quantitative studies of the CBS 
solution. 

IV. A QUANTITATIVE STUDY OF SHAPES IN THE CBS SOLUTION 

All the closed shapes (ii)-(vii) presented in last section can be grouped into two branches: The prolate ellipsoid 
including (ii)-(iv) and oblate ellipsoid branch including (ii), (v)-(vii), and these two branches are bifurcated from a 
sphere (ii). It is easily understandable if looking into the CBS solution (R) sin?/; = t/Rq + c$r lnr. When Rq > and 
Co = 0, this form gives nothing but sphere of radius Rq. Then when cq is a small quantity, positive cq leads to oblate 
ellipsoid and negative Co leads to prolate ellipsoid respectively. What if cq is large is not so evident and needs some 
reasoning. 

For analyzing how Co affects the shapes and the intervals of r such that | sin ip(r\ < 1, we use the parameterization 
sin-0 = r + cor lnr, and sketch the relation of simp with positive Co vs. r in FIG. 3. 

From this figure, we see that sinip = r + cor lnr is a single- valued and monotonous function of r; and it reaches its 
extremum at point r 2 which satisfies 

dsin ^ 1 , , i n • ( 1 + c 0\ /-,«•, 

— = 1 +c + c mr 2 = 0, i.e., r% = exp( ). (16) 

ar Co 

The extremum value of simp (r 2) is 

sin ip(r 2 ) = -r 2 c Q = -c exp( ), (17) 

Co 

which is negative for positive Co and vice verse. Taking the positive Co as an example in FIG. 3, corresponding to 
different Cq, shapes can appear in different intervals of r. There are three different cases of Cq € [0, 00). i) When c is 



critical such that sin-0(r 2 ) = — 1, we have c = c r rs 3.59112 from Eq.(17). The meridian principal curvature vanishes 



at point r 2 as c 2 = dsimp/dr = 0: it is the unphysically infinitely long capped cylinder with radius r 2 = l/c = 0.278 



from Eq. (17). ii) When Co < c r , we have simpfa) > —1 and shapes appear in one interval [0,7-4] in which r4 are 
the equatorial radii of the shapes such that simp(r^) = 1. Along the increase of r from to r4, the tangent angle 
ip(r) decreases from zero and reaches its minimum at point r 2 , then increases monotonously and reaches its maximum 
value tp — 7r/2 at point r4. These shapes certainly belong to the oblate ellipsoid branch, and their normals point 
outwards from Eq. (0). iii) When co > c r , we have sin7/>(r 2 ) < — 1, shapes can appear in two distinct intervals 
[0,ri], hi which r\ are the equatorial radii of the shapes such that sin^(ri) = — 1, and [7*3^4], because in both 
intervals | sin-0(r) | < 1 is satisfied. In interval [0, ri], the tangent angle ip(r) decreases from zero and terminates at 
ip(ri) = — 7r/2. These shapes certainly belong to the prolate ellipsoid branch, but the normals point inwards since 
— 1 < simp(r) < holds for whole interval r S [0, 7"i]. In distinct interval [7-3,^4], numerical studies show that only 
the self-intersection nodoidlike cylinder appears. 

In fact, using the parameterization simp = t/Rq + corlnr, each shape with a Co outside the domain (— l,c r ) can 
be found an identical shape with a value of Co within the domain (— l,c r ). For demonstrating this fact, we resort 
to the second form of the CBS solution sin tp = corln(r/r m ) in Eq. (R). Two parameterizations corln(r/r m ) and 
— cor ln(r/r m ) represent the same shape, but one's normal is opposite to another's from Eq.fl). If requiring that the 



normal points outwards, one of these two parameterizations must be eliminated. Furthermore we have a simplified 
form 

aiaip = c rln(r/r m ) 

= (c r m )(r/r m )ln(r/r m ) 

= c' xlnx, (18) 

where x = r/r m is the dimensionless length, and c — Cor m = coexp(—l/(coRo)) is a dimensionlcss spontaneous 
curvature. One can find that the c in two separate domains c € (—00, — e) and c G [0, e), in which e is the base of 
the natural logarithm, suffice to give all possible shapes with normals pointing outwards in the CBS solution. When 
c = ±e, the shapes are unphysically infinitely long capped cylinders. From the relation c — coexp(—l/(caRo)) and 
letting Rq — 1, all shapes can be mapped into a single domain cq G (—1, c r ) in parameterization sim/> = r/ Rq+cqt lnr. 
This is why there are only a few of spontaneous curvature Cq within a narrow domain in the last section are sufficient 
to give all possible types of shapes in the CBS solution, i.e., why each shape with a Co outside the domain (— l,c r ) 
can be found an identical shape with a value of cq within the domain (—1, c r ). 

In addition to the above qualitative analyses, numerical method will be used to characterize quantitatively the 
shapes in whole domain c$ G (—00,00). We introduce two ratios, semi-axis ratio z(ro)/ro and reduced radius ratio 
r v /r s . The so-called semi-axis z(r ) is defined by the half of the distance of two poles of a vesicle in the symmetrical 
axis z, which is 

pro 

z(r Q ) - z(0) = - / tan^(r)dr, (19) 

Jo 

where r$ is the equatorial radius of the vesicle satisfying smip(ro) = ±1. The ratios z(ro)/ro are plotted by dashed 
lines in FIG. 4. From FIG. 4(b), the cq approaches —1 and c r , the ratio z(ro)/ro approaches infinity: the shapes 
are infinitely long capped cylinders with radii 1 and 0.278 respectively. When Co increases from —1 to c r , the ratio 
decreases monotonously. We have in sequence the capped cylinder, the prolate ellipsoid and the sphere of unit radius 
at Co = 0, the oblate ellipsoid with a small positive Co, the CBS and the center of the CBS touches at Co = 2.4288, then 
self-intersecting inverted CBS. From FIG. 4(a) and FIG. 4(c), we clearly see that when cq < —1 or Co > c r , the ratios 
z(?"o)Ao are greater than 1. We have in sequence the capped cylinder, the prolate ellipsoid and the quasi-sphere with 
increasing |co|. When |co| approaches infinity, the constraint Isini/^r)! < 1 means that r\, cf. FIG. 3, must be very 
small. Thus we have approximately simp(r) ~ Tco?" for positive and negative cq respectively using L' Hospital method. 
It implys that we will have sphere again. All these shapes are closed and appear in intervals of r including the point 
r = 0, which are physically interesting and are plotted in the FIG. 4. When cq < —1 or cq > c r , in intervals [r^, > 0, r^] 
(cf. FIG. 3), we have the self- intersecting nodoidlike cylinders, which are less physically interesting and have not been 
plotted in the FIG. 4. As pointed out in |12), along with the increasing |co|, the number of self-intersections in unit 
length increases and it approaches a quasi-torus as |co| approaches 00. In fact, when —1 < cq < 0, the self-intersecting 
nodoidlike also appear in the intervals [r^ > 0, r^], but along with the decreasing |co|, the number of self-intersections 
in unit length increases. 

The second ratio is r v /r s , in which r v is defined by the radius of the sphere having the same volume V of a vesicle 
and r s the radius of the sphere having the same area A of the vesicle. The volume V and area A are 

J a cos tp 

f ro r 

A= 4tt / dr. (21) 

Jo C0S V 

Undoubtedly, we have r v /r s < 1 and the equality holds for sphere only. The ratio is plotted by dot-dashed lines in 
FIG.4. 

The scale invariant cor s is plotted by thin solid lines in FIG. 4. In FIG. 4(c), we plot both cor s (above the Cq 
axis) and — CQr s (below the cq axis ) by the same thin solid lines, in order to see the fact that each shape with a 
Co G (—00, —1) is identical to a shape with a Co G (c r , 00) and vice verse. But the normals of shapes parameterized 
by sin?/; = r + cqt lnr with Co G (c r , 00) point inwards. 

From both the qualitative and quantitative studies in this section, we can draw a conclusion that all possible closed 
shapes in the CBS solution (8) are, along with Co increasing from — 1 to c r , the capped cylinder, the prolate ellipsoid, 
the sphere, the oblate ellipsoid, the CBS, and the self-intersecting inverted CBS. 



V. THE SCALE INVARIANCE OF THE ENERGY AND THE SHAPE WITH MINIMUM ENERGY 

An important property of the local curvature energy is its scale invariance. This energy does not depend on the size 
of the vesicle but only on its shape [jlO| . If R is a solution with a local curvature energy, the rescaled shape R — ► R/fc 
with k > 0, and consequently (c\, C2, Co, dA) — > (feci, kc2, fcco, dA/k 2 ), is also a solution with the same local curvature 
energy. Two parameterizations of the CBS solution (8) having manifest scale invariance are 

sm.ip(r) = t/Rq + co?'ln(r/'ro), and sinip(r) = Corln(r/ro), (22) 

in which tq is a constant with length unit as r or Rq. It is evident that the two products c$r s and cqt v are two 
scale invariants independent of a particularly chosen parameterization. We used two particular parameterizations, 
simp(r) = r/Ro + corln?" and sin^(r) = corlnr, to search for the minimum energy shape. Both gives the same result 
cqt s = 1.04 and cor v — 1.00. The minimum energy is 0.480; hereafter we take the energy of sphere 8irk as the energy 
unit. Surprisingly, the CBS with touching center at cq = 2.4288 has energy 1.00. In FIG. 4 (b), the thick solid line 
shows that there is a energy minimum shape with cq = 1.200: it is a CBS. In general, the oblate ellipsoids including 
the CBS have lower energy than other shapes including the sphere, prolate ellipsoid and the self-intersecting inverted 
CBS have. 

VI. A COMPARISON OF OUR RESULTS WITH THE PREVIOUS EXPERIMENTAL AND 

THEORETICAL RESULTS 

Even all the previous numerical studies concentrated on the solution of equation when C = of Eq. (5), the obtained 
shapes |l(J,|l6[ included all the eight shapes above. And also, all these eight shapes except the self-intersecting cases 
have been observed in laboratory |16|-|l8|. 

Since our solution belongs to the situation C ^ of Eq.(5), our results coincide with the known theoretical 
results in one respect but differ from in the other. First, our approach supports the conclusion that the prolate 
ellipsoids may have higher energy than the oblate ones have Q. However, the standard instability analysis of a sphere 
starts from a slightly deformed sphere of parameterized form r = Rq + ^aj m Y] m (0, 0) where ai m is a set of small 
parameters corresponding to spherical harmonics Y/ m ($,</>). For small a; m , for example the only nonvanishing a2,o 
satisfying |a 2 ,o| << -Ro/2, the ellipsoids have C°° contour curves. But our ellipsoids have only C 1 contour curves, 
We must stress that the C 1 continuity suffices to ensure that the membranes are free from any force acting on the 
any point, otherwise the contour curve may be C . The second derivative of z(r) with respect to r is singular at 
point r = 0. One may feel uneasy of this singularity. In fact, it relates to energy density only, and the total energies 
(l/2)fc J(ci + C2 — co) 2 dA are limited as shown in FIG. 4. Second, the infinitely long capped cylinders appear when 
Co = —1 and Co = c r , and they correspond the infinite energy states, as shown by thick solid lines in FIG. 4 (a)-(c). 
From our results, the phases between the two sides of Co = —1 and Co = c r are separate. The latter c r distinguishes a 
oblate ellipsoid and a prolate ellipsoid phase; it is the same situation in the known phase diagram JlQfl . However when 
Co changes from Co < c r to Co > c r , the normal changes its direction from pointing outwards to pointing inwards. 
It seems that the membrane takes a global flip-flop procedure. But no such procedure takes place when Co changes 
from Co < —1 to Co > —1 even the transitions are also discontinuous. Third, there is an interesting number appeared 
both in our results and the previous ones, coi?o = 1-200. We see that from the FIG. 4(b) when cqRq — 1.200 where 
Rq = 1 in our parameterization, the CBS has the minimum energy. In the usual instability analysis of a sphere via 
infinitesimal deformation, the infinitesimally deformed oblate shape has the lower energy and more stable than the 
infinitesimally deformed prolate one whenever coi?o < —1-2 Q. 

VII. CONCLUSION 

An analytic solution for Helfrich spontaneous curvature membrane model [plpdj, which has a conspicuous feature 
of representing the circular biconcave shape, is systematically studied in this paper. Results show that the solution in 
fact describes a family of shapes, which can be classified as: i) the flat plane (trivial case), ii) the sphere, iii) the prolate 
ellipsoid (FIG. 2 (a)), iv) the capped cylinder (FIG. 2(b)), v) the oblate ellipsoid (FIG.2(c)), vi) the CBS (FIG.2(d)), 
vii) the self-intersecting inverted CBS (FIG. 2(e)) and viii) the self-intersecting nodoidlike cylinder (FIG. 2(f)). All 
these shapes have been found in numerical solutions of Eq. (5) with C — I. Except the self-intersecting cases, 



they all have real correspondence in vitro and in vivo on vesicle shapes 16 -lq]. The closed shapes ii)-vii) form two 



separate prolate and oblate ellipsoid branches which are bifurcated from a sphere. The oblate ellipsoids including CBS 
have lower energy than the prolate ellipsoids including sphere, and a CBS with c$r v = 1.00 has the lowest energy. 
The usual instability analysis of sphere leads to that the least stable shapes are prolate and oblate ellipsoid branches 
P,fl9[ , but it is the first time to give an explicit parameterization to show how these two branches come out gradually 
and analytically. 

We are indebted to Professors. Peng Huan-Wu and Zheng Wei-Mou for enlightening discussions. This subject is 
supported by National Natural Science Foundation of China. 
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FIG. 1. Sign convention. Four arrows mean positive directions for rotational axis z and radial axis r, tangent angle ip and 
arclength s, respectively. At the north pole both the arclength s and the tangent angle ip take zero values. 




(a) 



(b) 







(c) (d) (e) (f) 

FIG. 2. The nontrivial shapes in the CBS solution. The prolate ellipsoid (a), the capped cylinder(b), the oblate ellipsoid (c) 
the CBS (d), the self-intersecting inverted CBS (e) and the self-intersecting nodoidlike cylinder (f). 




FIG. 3. The relation of s'mip(r) vs. r and the intervals where the shapes appear. When < Co < Or, — 1 < smip{r2) 
shapes appear in [0,r4]. When Co > c r , sini/>(r2) < — 1, shapes appear in two distinct intervals [0, ri] and [r3,r4]. 
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FIG. 4. The energy (thick solid line), the scale invariant cor a (thin solid line), two ratios z(ro)/ro (dashed line) and r v /r s 
(dot-dashed line), in three domains — 7 < c < —1.2 (a), —0.98 < cq < 3.5 (b), and 3.6 < Co < 7 (c), for closed shapes. The 
bending elastic modulus k used to scale the energy of a sphere to unit. All quantities in four curves are dimensionless. In 
order to see clearly that fact that a shape with an exclusive value of Co G (— oo, —1) and a shape with an exclusive value of 
co G (c r ,oo) are identical except that the normals point inwards and outwards respectively, in last figure (c), we plot both 
cor a (above the Co axis) and — Cor s (below the Co axis ) by the same thin solid lines. We can also see that all shapes with 
co G ( — oo, — 1) and Co G (c r ,oo) have their corresponding identical shapes with Co G ( — 1,0). A single domain Co G ( — l,c r ) 
thus suffices to give all shapes. 
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